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' The gravitational collapse of a star is a warmly discussed but still puzzling problem, 
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problem, and reach the results: (I). The comoving coordinate system for the stellar 
system is only compatible with the zero-pressure free falling particles. (II). For the 



CJ ! free falling dust, there are three kind of solutions respectively corresponding to the 

oscillating, the critical and the open trajectories. The solution of Oppenheimer and 



Snyder is the critical case. (III). All solutions are exactly derived. There is a new 
kind singularity in the solution, but its origin is unclear. 
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^ ■ The gravitational collapse of a star is a warmly discussed but still puzzling problem, which 

H ; 

5r , not only involves the dynamics of the gases, but also the subtle coordinate transformation. 

The most enlighten model is the free falling dust model analyzed by J. Oppenheimer and 
H. Snyder [1]. In the paper, they said: "To investigate this question we will solve the field 
equations with the limiting form of the energy-momentum tensor in which the pressure is 
zero. When the pressure vanishes there are no static solutions to the field equations except 
when all components of TjJ' vanishes", "we are supposing that the relationships between the 
do not admit any stationary solutions, and therefore exclude this possibility." So they 
actually assumed the singularity exists, and ignored the dominant gravitational potential 
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itself is a more powerful source of pressure than the electromagnetic interaction [2, 3]. The 
gravity is conservative, so the free falling particles will move outward when they have passed 
across the stellar center. In addition, the global validity of the comoving coordinate system 
was employed without a careful check. Actually, the global simultaneity is not a general 
valid concept [4, 5]. 

In this letter, we give some more detailed investigation on this problem. From detailed 
calculation, we find that: 

1. The comoving coordinate system for the stellar system is only compatible with the 
zero-pressure free falling particles, so it may be not a globally valid coordinate system 
for the case. 

2. Similarly to the cosmological model, for the free falling dust we also have three kind of 
solutions, namely, the oscillating, the critical and the open trajectories. The solution 
of Oppenheimer and Snyder corresponds to the critical one. 

3. All exact solutions to the free falling dust can be manifestly derived, but the origin of 
a new kind singularity is not clarified. 

4. The relative equation between the normal spherical coordinate system and the co- 
moving one is simplified, but it is a nonlinear hydrodynamical equation, so the global 
solution might not exist in the general cases. 

The line element in the comoving coordinate system is given by 

ds^ = dt^ - e'^dr^ - e'^idO^ + sin^ 9d^^). (1) 

In this case, we have the 4-dimensional speed of the fluid = (1,0,0,0), and then the 
energy-momentum tensor becomes 

T^^ = diag(p, -Pe'^, -Pe"', -Pe"" sin^ 6). (2) 

In such coordinate system, we find all exact solutions can be solved by a bootstrapping 
method. That is, we firstly express all variables as functions of w and its derivatives, and 
then solve a simplified equation for w. 

By the continuity equation U^T^^^ = 0, manifestly If^d^p + (p + P)U'^.^ = 0, we get 

Idtw + dtw + = 0. (3) 
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For the Barotropic fluid, defining 
then we can solve w from (3) 

w = -2{w + B)+a{r), (5) 

where a is a function to be determined. Substituting (5) into the Einstein equation Gtr + 
hiTtr = 0, (k = SttG), namely [1, 6, 7] 

d^^w + ^{dtwdrW — dtwdrw) = 0, (6) 



we get 



B = -^w-\n{drw)+p{r), (7) 



where f3{r) is also a function to be determined. Substituting B into (5), we get w expressed 

by w 

w = w + 2ln{drw) - \n[A{l + -f)], (8) 

in which 7(r) = ^e^^"" — 1 is introduced for convenience. The effect of 7 is similar to that 
of the dimensionless energy density Q in cosmology as shown below, which determines the 
different evolving fate of the star. 

Substituting the above results into Gu + kTu = 0, we get p expressed by w 

p = [{-fdrw + 2dr-f)e~''" + dtw{3/4dtwdrw + dlw)] . (9) 

Taking dtp as variable independent of p and substituting (9) into Grr + i^T^r = 0, we get dtp 
expressed by w 

^tP = Tw^T — ^i'^^tr'^ + 'idtwdrw){dfwdrW - dtwdl^w + 2(9^76"""). (10) 
Substituting (9) into (10), we get the equation of consistence 

d%^w = -"^dtwdtrw + 9,(7e-"). (11) 
Integrating (11) with respect to r, we get 

d',w = -\{dtwf + ^e-^ + 5{t), (12) 
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where S{t) is a function to be determined. In Oppenheimer and Snyder's work, they assume 
7 = 5 = 0. As calculated below, their solution corresponds to the critical case. 
Substituting (11) and (12) into Gee + K,Tee = 0, we see the physical meaning of 6 

S{t) = -kP. (13) 

Relation (13) has an important implication, that is, the metric (1) is unsuitable to describe 
the realistic process of the star contraction, because in this case we definitely have drP 7^ 0, 
which contradicts (13). The comoving coordinate chart is usually a local coordinate chart in 
heavily curved space-time, and only in cosmology (1) is globally valid, because in this case 
P = P(t). This shows the limitation of the comoving metric (1). 

By (13), we find the only consistent matter in the space-time describe by (1) is the fluid 
with zero pressure. Of course, this model is extremely idealized, because the particles passing 
across the center will continue to move outward, and the collision among particles leads to 
high pressure and temperature. Such high pressure will balance the contraction[2, 3]. 

Now we solve (12) with 6 = —kP = 0. Setting w = |ln(-u), by (12) we get 

d^u = ^u-K {d,ur = l{s + jJ), (14) 

3 

where £{r) is a function to be determined. Again setting u = , i.e. making transformation 

w = ln(i;^), V = 62", (15) 

we get 

(d.vr = (16) 

Substituting (16) into (9), we learn the physical meaning of S 

drS = Kpv^dr-v, £{r) = £q + I Kpv^drvdr, (17) 

Jo 

where £q is a constant. £{r) corresponds to the comoving mass in the ball of radius r. 
Therefore, the quantity pv^dj-vdr is conserved. By their physical meanings, we certainly 
have drV > and dr£ > 0. From (17), we flnd v can be regarded as new spatial coordinate 
to replace r. Substituting drV into (8), we get 

= p= — (18) 
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Among the interim parameters {a, /3, 7, S, S), only (7, S) are independent functions. Their 
relations are given by 



;3 = ln^, (5 = lnp), a = 21n- 



6 = kP, 



n + 7 

where /3 is derived by substituting (16) into (7), and a is derived by the definition of 7. 
If 7 = 0, we get the solution derived in [1] as follows 



(19) 



v = ^i/l8S\t-to{r)\^. 



The curvature at t = to{r) is singular. Substituting (20) into (8), we get 

\drS{t-to) -2Sdrto\ 
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^12S'^\t-to\ 
Substituting (20) and (21) into (18), we find 



(20) 



(21) 



P 



(22) 



3K^12£^{t-toy 3K\[drS{t-to)-2Sdrto]{t-to)\' 

At the surface t = to{r) we have a singularity, which corresponds to the matter on this surface 
is passing through the stellar center. This space-time is bouncing. The phase trajectories 
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Figure 1: There are three kind of trajectories for the free falling dust. If 7 < 0, which means the 
kinetic energy of the particles is small, the orbits is oscillating. If 7 > 0, the orbits is open 

of (16) are displayed in Fig(l). As shown below, if 7 < 0, the space-time is oscillating, and 
if 7 > 0, the space-time is bouncing. 



6 



In the case 7 < 0, letting 7 = — ^, then we get the solution 



1 2 rf — 2v /—r^ r v^i 1 /r.n\ 

-?7 arccos \^ v[j]^ - v) = t - to • (23) 

The space-time is oscillating. Solving d^v from (23) and substituting it into (8), we get 

= V''"(^ ~^ ^=\\{ridrE - 6^9^r/)|t - tol - 2sgn(t - to)^^5rto] + ^4Ev^ridrr\\. (24) 
2r]V\/E{rf- — S) 



In the case 7 > 0, letting 7 = ^, we get the solution 



V 



V'ln— '-^== + ^v{r^^ + y)= \t-to\. (25) 

V'" + yV ^ 

The space-time is bouncing. Solving drV from (25) and substituting it into (8), we get 

= ^ ^' )=\[{r]drS - 6Sdrr])\t - to I - 2sgn(t - to)Sridrto] + 4v^wVr^l- (26) 

When t ^ to, the variables have the following asymptotic behavior in all cases, 

V = e"^^ ^^18£|t-to(r)P. (27) 

w \drS{t — to) — 2Sdrto\ , . 

6 2 , — . (28) 

■^'12^2|^-tolv/r+7 



^d..£VTT^ (29) 



3K\[drS{t-to) -28drto]{t-to)\' 

From the above calculation, we find the singularity at t = to implies the free falling 
particles pass across the stellar center, so it is a singularity caused by simplification of the 
matter model, but not a real one. The dynamics for the pressure-free fluid is ill-posed in 
mathematics, which results in the zero speed of sound and '5- wave solutions'. So we can not 
take such unrealistic model and solutions too serious, it just provides some intuitions. 

There are another kind singularity hides in the solutions. In the critical case 7 = 0, (21) 
and (22) become singular at time 

If to('") is not a constant, then t ^ to, the singularity is different from the one mentioned 
above, to(^) = const means a moving mass shell. So whether this singularity is caused by 
the globally invalidity of the comoving coordinate system or by the ill-posed model is not 
clear. Strangely enough, the function or f is normal at this time. 
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In the comoving coordinate system, the spatial coordinates are defined by the trajectories 
of a set of particles, and the temporal coordinate is defined by the proper time of the 
particles[l, 6, 7]. In cosmology, since the space is isotropic and homogeneous, this coordinate 
system works quite well. However in the heavily curved space-time, the Gaussian normal 
coordinate system is usually valid locally, and the global simultaneous hypersurface can not 
be arbitrarily defined[l]. (13) shows that, the comoving coordinate system may be only valid 
locally. In what follows, we examine the transformation between normal coordinate system 
and the comoving one. 

In general, the line element of an evolving space-time with spherical symmetry is given 

by 

rfs^ = bdt^ + 2fdtdr - adr^ - r^dO^ + sin^ Odif'^), (31) 

where (a, b, f) are smooth functions of (t, r). Making a transformation for t by solving a first 
order linear differential equation, we can remove / from the metric [6]. This is equivalent to 
set / = in (31). Since for a normal star, the metric functions (a, b, f) are at least continuous 
functions with bounded first order derivatives [3], then the solutions of the first order linear 
differential for the coordinate transformations, such as the one given in [6, 8, 9], have 
bounded second order derivatives, so they provide globally valid coordinate transformations. 
But the following analysis shows we usually can not find the global comoving coordinate 
system in an arbitrary regular space-time. 

For any smooth metric functions (a, b) with d^-b ^ 0, assume the following transformation 

T = T{t,r), R = R{t,r) (32) 

makes a comoving coordinate system with line element 

ds^ = dT^ - e^dR^ - e'^idO^ + sin^ 9d^^). (33) 

Substituting (32) into (33), and compare the results with (31) while / = 0, we get the 
equations for (T, R) as follows 

{{drTy + a)dtR - dr-TdtTdrR = 0, (34) 
OrTdtTdtR - {{dtTf - b)drR = 0. (35) 

They are a linear equation system for parameters {dtR,drR). The sufficient and necessary 
condition for nonzero solution is the determinant of the coefficient matrix vanishes, then we 
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have 

a{dtTy - b{drTf = ab. (36) 
Let u = —dtR/drR, substituting it into (35) and combining the equation with (36), we get 

= 0,T = (37) 



y/b — au^ \/b — 



au 



In (37), the signature ± can be removed by a transformation u —u, so we only need to 
consider the positive case. The integrable condition is given by 

dta 3 / dj-b dra\ 2 f dtb dta\ drb 



dtu + udr-u = -^u^ + ^ —]u'+[—, -]u-—. (38) 

2b \ b 2a J \2b a J 2a ^ ' 

Each solution u defines a comoving coordinate system. However, this is a 1+1 dimensional 
nonlinear hydrodynamical equation, which usually has only local smooth solution. The 
solution usually blows up in a finite time even for the smooth functions (a, b). 

In summary, from the above calculation and analysis, we find that the gravitational 
contraction is not so simple as it looks like. The comoving coordinate system for the stellar 
system is only compatible with the zero-pressure free-falling particles, so it is not globally 
valid coordinate chart. The valid coordinate transformation should be a 1-1 mapping with 
at least bounded second order derivatives in its domain of definition[8]. The free falling dust 
is just an idealized model, we can not treat it too seriously. The space-time itself becoming 
singular lags behind the matter source becoming singular, and then the conservative property 
of the gravitational potential may resist the formation of singularity [8]. 
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